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ABSTRACT 
The dark matter distribution in the very inner region of our Galaxy is still in debate. In the N-body 
simulations a cuspy dark matter halo density profile is favored. Several dissipative baryonic processes 
however are found to be able to significantly flatten dark matter distribution and a cored dark matter 
halo density profile is possible. The baryons dominate the gravitational potential in the inner Galaxy, 
hence a direct constrain on the abundance of the dark matter particles is rather challenging. Recently, 
a few groups have identified a tentative 130 GeV line signal in the Galactic center, which could be 
interpreted as the signal of the dark matter annihilation. With current 130 GeV line data and adopting 
the generalized Navarro-Frenk- White profile of the dark matter halo, for local dark matter density 
po = 0.4 GeV cm -3 and r s = 20 kpc we obtain a 95% confidence level lower (upper) limit on the inner 
slope of dark matter density distribution a = 1.06 (the cross section of dark matter annihilation into 



gamma-rays (av) > 



1.3 x 10 



17 cm 3 s 1 ). 



Such a slope is consistent with the results of some 



N-body simulations, and if the signal is due to dark matter, suggests that baryonic processes may be 

unimportant. 

Subject headings: Dark matter— Gamma Rays: general— galaxies: structure 



1. INTRODUCTION 

In the leading cold dark matter (CDM) model, struc- 
ture forms hierarchically bottom-up, with DM collapsing 
first into small halos, which then accrete normal matter, 
merge and eventually give rise to larger halos. Galaxies 
are thought to form out of gas which cools and collapses 
to the centers of these DM haloes (e.g., White & Rees 
1978). As shown in the high-resolution N-body simula- 
tions, the density profiles of CDM halos can be reason- 
ably well described by an universal form, independent 
of the halo mass, and the cuspy density profiles such as 
Navarro- Frenk- White ( NFW, Navarro et al. 1997) and 
Einasto (jEinasto 1119651 ) are found to be favored. The ef- 
fect of the baryons ignored in most previous simulations, 
however, are still to be figured out. One possibility is 
the so-called adiabatic contraction, i.e., when the Galaxy 
formed and the baryons contracted towards the centre, 
dark matter particles are pulled inward and their central 
density increases (e.g., Blumenthal et al. 1986; Gncdin 
et al. 2004). The resulting DM profile of the galaxies is 
expected to be more cuspy. Contrary to adiabatic con- 
traction, other baryonic processes, such as the gas bulk 
motions, possibly supernova-induced in regions of high 
star formation activity, and the subsequent energy loss 
of gas clouds due to dynamical friction, can transfer en- 
ergy to the central dark matter component or/and in- 
duce substantial gravitational potential fluctuations and 
finally give rise to a subsequent reduction in the central 
dark matter density (e.g., Navarro et al. 1996; Mo & 
Mao 2004; El-Zant et al. 2001; Mashchenko et al. 2006; 
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Ogiya & Mori 2011; Maccio et al. 2012; Pontzen & Gov- 
ernato 2012). The resulting DM profile of the Galaxies 
is likely cored. 

As shown in the recent high-resolution cosmologi- 
cal hydrodynamical simulations performed by Maccio et 
al. (2012; see Fig.l therein) and by Pontzen & Gov- 
ernato (2012; see Fig. 5 therein), with and without the 
effects of dissipative baryonic processes the inner dis- 
tribution of dark matter in Milky- Way-like objects are 
indeed rather different. Nevertheless these dark matter 
profiles only represent the mean of all simulated halos 
for a given mass at a given redshift. The scatter with 
respect to these mean values arises plausibly due to the 
different halo formation histories and due to the evo- 
lution of the expanding Universe (e.g., Navarro et al. 
2010; Wu et al. 2013). The average values may be un- 
able to describe accurately the DM halo of the Milky 
Way since its formation and evolution may not follow a 
prototypical spiral galaxy. In view of such uncertainties, 
observational data is highly needed to reliably constrain 
the inner structure of the Milky Way DM halo. How- 
ever, current microlensing and dynamical data can only 
rule out some extremely cuspy density profiles (Iocco et 
al. 2011) since in inner Galaxy, the gravitational poten- 
tial is dominated by the normal matter rather than the 
dark matter. It is well established that the prospect of 
detecting the products of dark matter particle annihila- 
tion is critically sensitive to the structure of the Milky 
Way DM halo. Moreover, a robust estimate of the cross 
section of dark matter particles annihilating into nor- 
mal matter (photons, electrons/positrons and so on) is 
not possible unless the DM halo profile has been reason- 
ably well determined. In turn, indirect DM gamma-ray 
searches can be used to study the distribution of DM in 
the innermost regions of the Milky Way halo as long as 
a positive signal (in particular an unambiguous gamma- 
ray line) has been identified. In this work, we adopt the 



hypothesis that the Galactic ~ 130 GeV line signal iden- 
tified in the publicly ava ilabl e Fermi-LAT dat a firstly by 
IBringman et al. 1 (|2012l ) and iWeniger I (|2012f) originated 
from the annihilation of dark matter. Under this hy- 
pothesis, we examine which dark matter density profiles 
are consistent with the signal, and then discuss the im- 
plication. 

2. CONSTRAINING THE DARK MATTER DENSITY 



PROFILE AND (ctV)j, 



WITH THE 130 GEV LINE 



SIGNAL 



2.1. The 130 GeV line signal 

High energy 7-ray line is of greatest interest in look- 
ing for the signal of dark matter annihilation. After 
analyzing the publicly avail able Ferm i -LAT 7-ray data, 
IBringman et al. I (|2012l) and IWeniger 1 (J2012D found pos- 
sible evidence for a monochromatic 7-ray line with en- 
ergy ~ 130 GeV. Later independent analyses confirmed 
such an excess. The center of the most prominent signal 
region is a round the Galactic centre Sgr A* with an of f- 
sct > 1.2° (jTempel et al. |[20ll[Su~&: Finkbeiner 1120121 ), 
which has been thought to be at odds with the DM ori- 
gin. However, within the DM annihilation scenario such 
an offset can be interpreted by the limited statistics of 
the c urrent 7— ray line s ignal consisting of only ~ 14 pho- 
tons (jYang et al. II2012J ). Furthermore, it is interesting 
to note that there is a small wiggle in the electron spec- 
trum of PAMELA/Fermi-LAT at energies ~ 100 GeV, 
which may be interpreted as being consistent with the 
130 GeV line signal (Feng et al. 2013). 

With a typical cuspy DM densit y profile such 

as Navarro -Frenk-White (NFW, iNavarro et al. I 

(|1997f )) and lEinasto I (|1965f ). annihilation cross sec- 



tion {<JV)y 



2 — 5 x 10 27 cm 3 s 1 is needed to 



produ ce the signal data (|Weniger 1 120121 : 1 Temp el et al. I 
I2012D . which however is larger than the upper limit 
(~ 10~ 27 cm 3 s _1 ) set by the non-detection of 130 
GeV line in other regions by a factor of quite a few. 
This puzzle has been taken to be a piece of evidence 
against the dark matter origin of the 130 GeV line in 
the Milky Way center (e.g., Huang et al. 2012). As 
already mentioned in Sec. 1, in the Galactic center the 
DM density profile is rather uncertain, it is thus not 
confidential to refute the dark matter model just based 
on the tension between the values or upper-limits of 
(cri>) xx _j. 77 inferred in different regions. 

In principle, one may be able to find an ideal region 
to reliably evaluate (crv) xx _» 77 . In such a region the fol- 
lowing conditions should be met, including (i) The dark 
matter density distribution is reasonably determined by 
astrophysical observations; (ii) The signal-to-noise is rel- 
atively high; (iii) The contribution of the DM substruc- 
ture is expected to be not dominated. For the Milky Way, 
the dynamical data plays a role in constraining the dark 
matter distribution for r > 3 kpc, where r is the distance 
to the Galactic center (e.g., Sofue 2012). Observation- 
ally there is still no strong evidence for the existence of 
abundant substructures of the Milky Way dark matter 
halo. The N-body simulation Aquarius found out that 
abundant substructures of the Galaxies might present at 
r > 20 kpc (Springcl ct al. 2008). The contribution 
to the J- factor may be non-ignorable (i.e., the contribu- 
tion is about the same as the smooth halo) any longer at 
ip w 30°. Unless for the very cuspy dark matter halo, the 



signa l-to-noise is also accept able in the region ip > 30° 
(e.g.. IBringman et al. 1 12012). In view of these facts, the 
region ip £ (20°, 40°) excluding b < 10° may be a suit- 
able region to constrain/measure (ctw) xx ^. 77 , where b is 
the Galactic latitude and ip is the angle to the Galactic 
centre. As an unbiased constraint on the dark matter 
density profile, the data in other regions should be taken 
into account, too. It should also be mentioned that 
the dwarf galaxies are also ideal candidates to constrain 
the annihilation cross section of dark matter. Geringer- 
Sameth & Koushiappas (2012) performed a joint analysis 
of dwarf galaxy data and found that the upper limit on 
the annihilation cross section to a two-photon final state 
is 3. 9^3.7 x 10~ 26 cm 3 s" 1 at 130 GeV (see also Huang et 
al. 2012), which is well above that needed to account for 
the signal identified in the Galactic center (see below). 

2.2. Fermi LAT data Analysis 

In last subsection we have suggested that ip £ 
(20°, 40°) excluding b < 10° may be a suitable re- 
gion to constrain/measure (crw) xx ^. 77 . In reality, for a 
given dark matter density profile the intrinsic (crv) xx ^ 77 
should be smaller than the values or upper limits inferred 
from any other regions. Currently the tentative 7— ray 
signal is present in a very compact region, hence we are 



only able to set an upper limit on (o~v) 



xx-m 



and 



pose a constraint on the slope of the dark matter density 
profile. 

For such a purpose, we analyze the publicly available 
Fermi-LAT data in several regions: including (a) ip < 2°, 
w hich covers the most p rominent signal region identified 
in lTempel et~ai~l ((20121 ): (b) the region iP £ (2°, 6°); (c) 
the region ?/> £ (6°, 10°);(d) the region ip £ (10°, 20°) ex- 
cluding b < 10°; (e) the region ip £ (20°, 30°) excluding 
b < 10°;(f) the region ip £ (30°, 45°) excluding b < 10°; 
(g) the rest region ip £ (45°, 180°) excluding b < 10°. 
We take into account the data in the time interval from 
4 August 2008 to 18 April 2012 (MET 239557417 - MET 
356439845), with energies between 20 and 200 GeV. We 
used the standard LAT analysis software (v9r27pl jj. To 
reduce the effect of the Earth albedo background, time 
intervals when the Earth was appreciably in field-of-view 
(FoV) , specifically when the center of the FoV was more 
than 52° from zenith, as well as time intervals when parts 
of the ROI were observed at zenith angles > 100° were 
also excluded from the analysis. The spectral analysis 
was performed based on the P7v6 version of post-launch 
instrument response functions (IRFs). ULTRACLEAN 
dataset was selected to avoid the contamination from the 
charged particle. The spectra are shown in Fig.l. 

To investigate the line signal at 130 GeV we use the 
unbinncd analysis met hod which is similar to the one de- 
scribed in (Ac kermann et aLl 120121 ). It should be noted 
that in drawing fig.l we binned the counts to different 
energy bins, which may introduce fake signal when the 
counts number is small. However, in fitting the data to 
derive the possible line signal or upper limits we use the 
unbinncd analysis, say, the likelihood function was built 
by multiply the probability distribution function of each 
photons in the assuming model. This method can mini- 
mize the fake signal due to binning and take advantage of 

4 http://fermi.gsfc.nasa.gov/ssc 
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Fig. 1. — The 20-200 GeV spectra of our Galaxy in seven regions 
described in the text (Please note that b < 10° is excluded for the 
last four regions). J 7 is the differential gamma ray flux. 

the full information of the observed data. The likelihood 
is described as 



C = UfS(E i ) + (l-f)B(E i ), 



(1) 



where S(Ei) and B{Ei) represent the signal and back- 
ground function, respectively, both are normalized to 1, 
and i runs over all the photons; / is the signal fraction 
and has been set to be in the range [—1,1] for line signal 
search and [0, 1] for getting upper limits; B(Ei) takes the 
form 

B(E. t ) ~ E- r e(Ei), (2) 

where e(Ei) is the exposure generated by the gtexpcube2 
routine. S{Ei) is derived by convolving the energy dis- 
persion function and exposure. Pyminuit are used to 
find the maximum of the likelihood. The MINOS asym- 
metric error at the level Aln£ = 1.35 is adopted to 
get upper limit corresponding to a coverage probability 
of 95%. To see the possible systematics due to the dif- 
ferent choice of the background spectrum, we redo the 
analysis by using two different background spectral tem- 
plate B(Ei): (1) adding an exponential cutoff to the pure 
power law spectrum and leaving the cutoff energy to be 
free, say, 

B(E t ) ~ E^exp^E./E^ME,), (3) 

(2) a log-parabola spectrum 



B(Bk) ~ ( 



Ej, 

E b ' 



-(a+/3log(E,/E b )) 



e(Ei). 



(4) 



We do not find significant improvement of the fitting, 
thus for simplicity we use the pure power law spectrum 
as our fiducial background spectral model. 

5 http://code.google.eom/p/pyminuit. 
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Fig. 2. — The energy dispersion used in our analysis. The form 
of energy dispersion is described in the text and has already been 
normalized. 

The region we choose is all larger than several de- 
grees, meanwhile, the angular resolution above 20 GeV is 
within 0.2 degree, which is much smaller than the size of 
the region of interest. Thus we neglect the point disper- 
sion function(PSF) in the analysis. On the other hand, 
the energy dispersion is extremely important in the line 
searching process. In this work we focus on the tenta- 
tive 130 GeV line, and adopt the energy dispersion at 
130GeV with the form described on the website 0. The 
P7ULTRACLEANV6 version of the instrument response 
functions(IRFs) is used and the final energy dispersion 
is averaged for different incidence angles. The derived 
energy dispersion used in the analysis is shown in Figj5] 
and the signal flux (or upper limits) obtained in different 
regions are summarized in Tab.l. 



TABLE 1 

The signal flux as well as the 95% confidence level upper 

limits obtained in different regions. please note that 

b < 10° is excluded for the last four regions. 



Region 


Position 


Signal/upper limit (fO 8 cm 2 s 1 sr v ) 


a 


ip < 2° 


3.8 ± 1.6 


b 


»/>6 (2°, 6°) 


0.55 


c 


4> G (6°, 10°) 


0.34 


d 


ip 6 (10°, 20°) 


0.133 


e 


V> e (20°, 30°) 


0.028 


f 


ijj e (30°, 45°) 


0.0084 


g 


i/> > 45° 


0.0015 



2.3. Constraints on (cru) xx _>. 77 and the dark matter 
density profile 

To constrain dark matter distribution in inner 
Galaxy, we implement spherically symmetric generalized 
NFW profiles 



PdmM = Ps{r/r s ) a (l + r/r s ) 



-3 + Q 



(5) 



where r s restricted in the range 10 — 35 kpc is the scale 
radius (jlocco et al. 1 1201 11 ) and a is the inner slope for 
the NFW profile (it is found in the N-body simulations 
that 0.9 < a < 1.2, however the baryon compression 

6 http://fermi.gsfc.nasa.gov/ssc/data/analysis/documcntation 
/Cicerone/Cicerone_LAT_IRFs/IRF_E_dispersion.html 
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Fig. 3. — The annihilation line flux expected in the theoretical 
models vs. the current 130 GeV line data. See the text for details. 

may give rise to a ~ 1.7). One main reason for adopting 
the generalized NFW profiles is that they can reproduce 
the three types of dark halo profiles presented in Maccio 
et al. (2012). The normalization of the DM profile is 
set by the local DM density, i.e., p = P dm (Rq) and 
p s = po(r Q /r s ) a (l + r Q /r s ) 3_a , where R Q « 8.5 kpc is 
the distance from the Sun to the Galactic center. The 
fiducial interval po = 0.4 ±0.1 GeV cm -3 is adopted in 
line with recent astrophysical measurements (Salucci et 
al. 2010). The 7-ray flux produced by DM annihilation 
can be written as 



$(Afi,£L 



1 
4- 



(ot;) xx _, 77 diV 7 
2m x d£ 7 



x J(Afi)Afi, (6) 



where m x is the mass of DM particles and diV 7 /d.E 7 



26 {E^ 



( ) is the differential energy spectrum of 7-rays. 



The astrophysical factor (J) is defined as 



J(Afi) 



1 
Aft 



dfi / dZp 2 (r(/)), 

AS! J LOS 



(7) 



where r(l) is the distance to the center of the object 
which is a function of line of sight (LOS) distance I. 

For region (a), there is a tentative gamma-ray line 
signal. If interpreted as dark matter particles annihi- 
lating into a pair of photons, the cross section of dark 
matter annihilation ((crw) xx _ >77 a ) can be inferred with 
equation ©. For regions (b), (c), (d), (e), (f) and (g) 
we have the 95% confidence level upper limits of the line 
flux and then get the constraints on the annihilation cross 
section. As long as 



Wxx-m,a < min{(o"u) xx _, 77ib , (crv) xx ^ 

\ CrU /XX->774' \ (7U /XX-!-77>oj 
\ aV )xX^llJ ' \ CrW /XX->77,g)' 



>77,g/' \°) 

the Dark Matter profile is in agreement with the 130 
GeV 7— ray line data. Such a fit to the current tentative 
130 GeV 7— ray line data suggests that a > 1.17 and 
(av) xx -^ 11 < 7.5 x 10~ 28 cm 3 s" 1 for r s = 20 kpc and 
po = 0.4 GeV cm" 3 (see the dotted line in FigEJl. 

To get more robust constraint, we adopt a com- 
bined likelihood analysis in all seven regions to constrain 
(cru) xx _ ! . 77 and a for r s ranging from 10 kpc to 35 kpc. 
The method is similar to that described in iTsai et all 



r s (kpc) 

Fig. 4. — The 95% confidence level lower (upper) limit on the 
generalized NFW distribution parameter a (((Tt;) xx _ >77 ) set by 
current 130 GeV line data for po = 0.4 GeV cm -3 . Following 
llocco et al. 1 1 120111 ). r B is restricted within the range 10 — 35 kpc. 



(J2013I) and the combined likelihood is calculated as 

L c = HLi, (9) 

where Li is the unbinned likelihood in the i— th region. 
The definition of the unbinned likelihood is similar to 
that in last section, but now the influence of (<r?;) xx _>. 77 
as well as a should be taken into account. For such 
a purpose we modify the signal ratio / in Eq.(l) to 
/ = $((CTw) xx ^ 77 ,a)/$ obs , where $((av) xx ^gj, a) is 
the flux predicted in the DM model (i.e., eq.(|6])) with 
varying (erv) xx _j. 77 and a, and $ b s is the observed 
gamma ray flux (i.e., the integration of J 7 in Fig.l in 
the whole energy range). In o ur approach the p rofile 
likelihood technique is adopted (jRolke et al. Il2005h . To 
constrain (cr^) xx _> 77 we treat a as a nuisance parameter 
and vice versa. Aln(L c ) = 1.35 is adopted to get up- 
per (lower) limit corresponding to a coverage probability 
of 95%. In this work we do not constrain po (or alter- 
natively p s ) since it couples with the annihilation cross 
section, i.e., (av) xx ^ 11 a 1/Po- Cur general 95% con- 
fidence level lower (upper) limit on a ((crw) xx _ >77 ) as a 
function of r s is presented in Figf4] For r s = 20 kpc and 
po = 0.4 GeV cm" 3 , the 95% confidence level constraints 
are a > 1.06 and (ctu) xx _ >77 < 1.3 x 10" 27 cm 3 s" 1 



respectively. The required (av) 



xx-m 



is consistent 



with the constraints set by the non-detection of a reli- 
able signal in dw arf galaxies, di f fuse Galactic halo and 
Gala xy Clusters (|Strigari 1 120121 : iBringman fc Weniger I 
l2012f l. E ven for r s = 10 kp c, the smallest value sug- 
gested in llocco et al. I (|2011f ). a > 0.96 is needed and 
hence the HFR dark matter profile suggested in Maccio 
et al. (2012) is disfavored, implying that the baryonic 
processes that can considerably flatten the central dark 
matter distribution might not play an important role in 
the evolution of the Milky Way. 

In Eris, a very recent high-resolution cosmologi- 
cal hydrodynamics simulation of a realistic Milky- Way- 
analog disk galaxy, the dark matter density profile within 
~ 1 kpc is found to be very flat and the peak of the dark 
matter density profile is ~ a few 100 pc away from the 
Galactic centre (Kuhlen et al. 2013). It is interesting 
to check whether the specific cored dark matter distri- 
bution found in Eris is consistent with the current line 
data or not. The fits to the data are shown in Figj3l 
The dashed line represents the best fit. In region (b) 



the divergency between the model and the data is by a 
factor of > 4. The dash-dotted line (i.e., the so-called 
Erisumit model) represents the case that the predicted 
flux of Eris model in region (a) has been normalized to 
a flux w 1.17 x 1CT 8 cm" 2 s~ 1 sr" 1 , i.e., the 95% confi- 
dence level lower limit of the line signal. In such a fit, 
the model predicted line flux in region (b) is still above 
the upper limit. We thus conclude that the Eris model 
has some tension with the line data. 

3. SUMMARY 

In the very inner region of the Galaxy, the dark 
matter abundance is much less than the normal matter, 
hence current micro-lensing and Galactic rotation curve 
data can not directly constrain the dark matter density 
distribution. The numerical simulations in principle can 
solve such a problem. However, our current knowledge 
of the history of the Galaxy is very limited and the com- 
plicated physical processes being able to shape the dark 
matter distribution are hard to fully address. That's why 
so far the distribution of dark matter in the very inner 
region of the Galaxy is still in heavy debate (e.g., Gncdin 
et al. 2004; Maccio et al. 2012; Pontzen & Governato 
2012; Kuhlen et al. 2013). Recently several groups have 
identified tentative 130 GcV 7— ray line which might be 
due to the annihilation of dark matter particles in the 
inner Galaxy. In this work we adopt the hypothesis 
that these signals are due to dark matter annihilation 
and use this hypothesis to examine which DM profile is 
consistent with the line. Our finding is that at the 95% 
confidence level, the dark matter density profile towards 



the center should be not shallower than r -106 (for the 
generalized NFW profile with r s = 20 kpc) and the dark 
matter annihilation cross section should be smaller than 
(cru) xx ^ 77 = 1.3 x 10~ 27 cm 3 s^ipo/OA GcV cm" 3 )- 2 
(see Fig. 4). Such a density profile is in agreement with 
that found in some N-body simulations and the baryon 
compression effect might play a role. The dissipative 
baryonic processes that are able to considerably flatten 
dark matter profiles seems to not play important roles in 
our Galaxy, implying that the formation and evolution of 
Milky Way may not follow a prototypical spiral galaxy. 
Finally we'd like to caution that the tentative 130 GeV 
7— ray line has not been officially confirmed by the Fermi 
collaboration yet. Whether our constraints on the dark 
matter distribution in very inner Galaxy and on the cor- 
responding annihilation cross section are robust or not 
will be directly tested by the upcoming pass 8 data of 
Fermi LAT, in which the amount of usual data at ener- 
gies greate r than 10 GeV is ex pected to be boosted by 
some 60% (|Bloom et al. I [20121) . 
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